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Intense laser ionization expands Einsteins photoelectric effect rules giving a wealth of phenomena 
widely studied over the last decades. In all cases, so far, photons were assumed to carry one unit of 
angular momentum. However it is now clear that photons can possess extra angular momentum, the 
orbital angular momentum (OAM), related to their spatial profile. We show a complete description 
of photoionization by OAM photons, including new selection rules involving more than one unit 
of angular momentum. We explore theoretically the interaction of a single electron atom located 
at the center of an intense ultraviolet beam bearing OAM, envisaging new scenarios for quantum 
optics. 



PACS numbers: 03.67.Mn, 42.50.Dv, 42.65.Lm 

During the history of Physics, Light-Matter interac- 
tion has been the fundamental path towards understand- 
ing new phenomena and testing some essential theories, 
as in the case of photoionization. Approximately until 
1992, all physical theories have described light according 
to three features: energy, linear momentum and polar- 
ization. The latter, which is purely related to the elec- 
tric field direction of the propagating light, yields an ef- 
fective light angular momentum. But, as recognized by 
Allen and coworkers [1], light possess another degree of 
freedom: the orbital angular momentum (OAM), which 
rather than being associated with polarization, it is re- 
lated to the spatial profile of light. This newborn de- 
gree of freedom has kindled a huge activity in different 
lines of research, ranging from micro and nanoparticle 
trapping [2HI] to quantum state engineering in Bose- 
Einstein condensates [5], multiphoton entanglement [6j- 
[5] for quantum information applications and molecular 
spectroscopy [5]. Recently, both femtosecond and high- 
power OAM beams have been generated experimentally 
using holographic plates [101 [IT] . In this work, we revisit 
the Einstein photoionization scenario [12] . but now tak- 
ing into account the orbital angular momentum of light. 
This allows us to unveil new phenomena beyond the stan- 
dard photoionization. 

Photoionization has attracted a broad interest both 
from a fundamental theoretical point of view and from 
the standpoint of applications. Laser photoionization, 
particularly strong field photoionization, has been a very 
active research topic over the last decades. Many new 
effects have been reported, such as ATI (Above Thresh- 
old Ionization), tunnel ionization, high-order harmonic 
generation, etc. P31 [14] Besides, much effort has been 
devoted in the recent years to the possible inhibition of 
photoionization at high frequencies and for very strong 
fields. Experimental activity on ultra-strong field ion- 
ization has recently reached the relativistic domain [15] . 
All this photoionization literature can be classified into 



two regimes: the electric-dipole regime (where the mag- 
netic field of the light can be neglected) and the non- 
dipole regime. In the electric-dipole regime, light beams 
carry the standard angular momentum (one h unit), and 
the selection rules avoid the possibility of one-photon ex- 
citation of atomic transitions with angular momentum 
variation larger than one H. However, one can over- 
come this limitation by considering multi-photon effects 
with very intense lasers. In the non-dipole regime, the 
light-atom interaction is more complex, exciting not only 
atomic transitions with angular momentum change equal 
to one unit H, but also atomic transitions with larger an- 
gular momentum variation. In this work we present for 
the first time photoionization with light beams carrying 
OAM which give rise to new selection rules out of both 
the electric-dipole and the non-dipole regime, opening 
new perspectives for atomic transition excitations. 

Much numerical work has been devoted to 3D pho- 
toionization studies. Strong field photoionization implies 
ionized electrons that escape at high energies, thus re- 
quiring large and dense numerical networks (to describe 
both electrons far away from the atomic core and en- 
ergetic electrons driven by the laser). When a linearly 
polarized laser light is used, numerical solutions are rela- 
tively simple due to the cylindrical symmetry of the prob- 
lem and only two dimensional numerical grids are needed. 
However, there exist relevant examples where such cylin- 
drical symmetry is not possible [16!, and a true three 
dimensional numerical grating is needed. This demands 
a much larger complexity of the numerical simulation. 
Here, we present an accurate description of the atomic 
photoionization induced by a beam carrying OAM. This 
scenario requires a true 3D simulation to fully account for 
both the three-spatial dimensions of the electron quan- 
tum state and the 3D-spatial profile of the beam (includ- 
ing the transverse profile related to the OAM). 

In this paper we address the interaction of a pulse 
beam carrying OAM with the simplest atom: hydro- 
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gen. This scenario provides the paradigm for fundamen- 
tal questions: Is the OAM of light transferred to the 
electron quantum state? There are still open questions 
about the angular momentum transfer between matter 
and light [2J [T7H2U] . The angular momentum of light can 
be separated in OAM and spin momentum in the parax- 
ial regime, but how are the OAM and the spin momen- 
tum of light transferred to the matter? In some previous 
works [T8H2"U] . where an ensemble of atoms was consid- 
ered, it was shown how the OAM can be coupled to the 
center of mass of the atomic ensemble. The same inter- 
action description out of the paraxial regime [19] is more 
plentiful. In contrast with previous works, we are consid- 
ering the interaction of a Laguerre-Gaussian intense laser 
beam with a single quantum state placed near its optical 
vortex, taking into account the general form of the quan- 
tum electromagnetic-interaction-field hamiltonian, with- 
out neglecting any magnetic term and without restricting 
to any multipolar approximation in the transverse plane, 
going beyond into the comprehension of the OAM light 
coupling with a quantum system. To clarify the picture, 
we investigate photoionization with OAM beams in the 
Schrodinger regime, using both analytical and numerical 
tools. 
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FIG. 1: Light-Matter Interaction Scheme. The ad- 
dressed problem consists of: A temporal pulse with a well- 
defined polarization and a transverse profile that takes into ac- 
count the OAM. The initial state of the electron corresponds 
to the fundamental state of the hydrogen atom. 
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where is the step function, a the Bohr radius, c the 
speed of light, lo the carrier wave frequency, and A Q 
the amplitude of the wave (it includes the polarization 
state). The transverse spatial structure of the pulse 
beam is accounted by the functions LGg iP (p, (f>, z;u>/c); 
the Laguerre-Gaussian modes pQ. They are character- 
ized by a width wq at z = (the beam waist), and by 
the indices I = 0,±1, ±2, ... and p = 0,1,2,..., rep- 
resenting the, so-called, winding (or topological charge) 
and the number of nonaxial radial nodes of the mode, 
respectively. Laguerre-Gaussian modes contain an az- 
imuthal phase e lt ^ which gives rise to a discrete OAM 
of £h units per photon along their propagation direction. 
The complete spatial spectrum of scalar wave fields pre- 
pared in arbitrary superpositions of Laguerre-Gaussian 
(or other paraxial) modes can be measured by using a 
simple interferometric scheme [21]. The fact that the 
associated electric field amplitude now depends on the 
transverse position, in contrast with plane waves, will be 
shown below to give rise to unexpected phenomena. 

We assume that the hydrogen nucleus is unaffected 
by the electromagnetic field. If so, the electron-field 
coupling evolution will be described by the following 
Schrodinger equation 
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where ip(r,t) is the electron quantum state, V(r) the 
Coulomb potential originated by the hydrogen nucleus, 
m the electron mass, q the electron charge, A(r,i) the 
vector potential, which in our case is given by expres- 
sion (JlJ, and p = — ifiV the linear momentum operator, 
satisfying the canonical commutation relation [x, p] = ih. 



I. LIGHT-MATTER INTERACTION SCHEME 

We begin by considering a pulse beam propagating 
along the z-direction with a temporal envelope wave pa- 
rameterized by a quadratic sinus (see Fig. [I]). This tem- 
poral envelope has a frequency uj e — ir/N cyc T, where N cyc 
and r are the cycle number and the period of the carrier 
wave, respectively. A hydrogen atom is assumed to be 
localized at the origin of the reference system and expe- 
riences a vector potential, associated to the pulse, of the 



II. SELECTION RULES WITH OAM 

We first extract a new set of selection rules for the in- 
teraction of atoms with OAM beams and point out the 
essential differences with pulse beams consisting of plane 
waves. The Hamiltonian (|2j) of the system can be ex- 
pressed in the usual form; H = Hq + Hz + Uxxi where 
Ho is the free part, whereas Hi = — q (p • A(r,£) + 
A(r,i) • p)/2m and Hxx = q 2 A 2 (r,t)/2m refer to the 
interaction parts. Representing the quantum state in a 
spherical basis ^(r) = J^l m u L,M{r)Yjf (6, <p) (all radial 
dependence is in the functions UL,Af(^), while the angular 
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dependence remains in the spherical harmonic functions 
Yjf{6,(p) instead), the first interaction contribution can 
be written as 

^f\M^i)=i\i.E i -E f )^ f \v-k{v,t)\^ i ) . (3) 

Here, Ei and Ef are the unperturbed energies of the ini- 
tial and final states, respectively. Taking into account 
the vector potential given by ([!]), and within the dipolar 
(A ^ a a ) and the transverse spatial (too 3> A) approxi- 
mations, we derive the following set of selection rules for 
beams carrying any arbitrary £ units of OAM (see the 
Appendix) : 

\AL\ < \£\ + 1, AM = £ ± 1, with AL + \£\ + 1 even, (4) 

where we have assumed that the quantization axis is 
along the beam propagation direction. In contrast with 
plane waves (for which £ = and |AL| = 1), signifi- 
cant variations of the angular momentum are to be ex- 
pected. In terms of photons, the selection rules Q can 
be conceived as the absorption of photons carrying a to- 
tal angular momentum j = £ + s in the propagation di- 
rection, where s indicates the polarization part (or spin 
momentum, s = ±1 for right- and left-circular polariza- 
tion) . We would like to remark that these selection rules 
originate from the transverse profile, despite the dipolar 
approximation. Moreover, the second contribution for 
the interaction Hamiltonian Hzx yields, in the case of 
plane waves, a constant term, producing a ponderomo- 
tive force [22]. In our case, due to the transverse profile 
of the beam, this Hamiltonian part produces two contri- 
butions. One acting as a ponderomotive force, while the 
other one, remarkably, gives rise to new selection rules 
(see the Appendix): 

\AL\ < 2\£\, AM = 2£, with AL even. (5) 

We should remark that the domain of applicability of 
selection rules Q and Q extends beyond the photoion- 
ization problem. 

III. HYDROGEN SIMULATIONS 

Selection rules Q and ^ constitute our first main 
result. To gain a more complete picture, we now pro- 
ceed with the exact description of the hydrogen electron 
state ionization by beams carrying OAM. Prior to the 
interaction, we assume the electron to be in the ground 
state, ipi(r) = 1/ ^J-Ka i e~ r l ao . We simulate the evolu- 
tion of the electron quantum state when the incoming 
pulse has N cyc = 3, £ = 1, p = 0, an angular frequency 
u) = 1 au (atomic units where K = m = q = 1, and 
ui = 2tt x 6.57 • 10 15 s _1 , ultraviolet), a period r = 2n au 
(152 as), and two possible polarizations: linear (in the 
x-direction) and left-handed. We choose a beam waist to 
satisfy the paraxial regime, w = 9 • I0 4 au (4.79 //m), 
which is much larger than the characteristic size of the 
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FIG. 2: Carrier Envelope Phase for a Laguerre- 
Gaussian Beam. In the main figure, the polarization of 
the pulse beam is represented for a transverse plane at z = 0, 
just when the pulse is maximum. The atom is centered at the 
origin (coincident with the beam vortex), where the electric 
field is zero. The arrow size is proportional to the electric 
field amplitude strength, which increases linearly in the ra- 
dial direction up to amplitudes of about 5 au in the boundary 
distance, 20 au (lnm), and it is harmonically modulated by 
the azimuthal position. The blue line is the nodal line of the 
electric field, and rotates with the carrier wave frequency w, 
setting the electric field distribution during the pulse inter- 
action. The temporal dependence of the pulse beam is also 
shown at four different positions of the transverse plane. No- 
tice that the carrier envelope phase (CEP) differs from the 
azimuthal dependence, not the radial. 



atom (wq 3> a Q ). Our atom is centered at the beam 
vortex, interacting with its vicinity, where the electric 
field amplitude is much weaker than the maximum one 
(reached at a distance Wq/V^)- This imposes the need 
for very intense lasers; in the proximity of the vortex 
singularity the electric field amplitude increases linearly 
(when \£\ — 1). For example, an electric field E ~ A cu 
of about 10 4 au (10 13 V/cm), would give rise to 5 au 
amplitudes (during the pulse peak) at distances of about 
20 au (1 nm) from the singularity. In order to clarify 
the structure of the electric field, in Fig. [2] we repre- 
sent the polarization of the electric field in the transverse 
plane z = 0, when the pulse achieves its maximum value. 
Figure [2] also plots the pulse beam with respect to time 
at four different positions in the transverse plane. No- 
tice the variation of the carrier envelope phase (CEP) 
depending on the azimuthal position. In fact, all the 
possible CEPs are encompassed in a circle around the 
singularity, in contrast with standard few cycle pulses 
technology [23] , where much effort has been done to lock 
the carrier-envelope offset. 

Since the electric field has a frequency w = 1 au (larger 
than the bound hydrogen energy 0.5 au), ionization is to 
be expected. In fact, after the first pulse cycle we verify 
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FIG. 3: Initial Ionization of the Electron Quantum 
State, (a) In the upper row, projection of the excited state 
onto the plane xy (J dz\Sip(x,y, z)\ 2 ) at four different times 
when the pulse beam is linearly polarized. The electron is be- 
ginning to be ionized in the first cycle of the pulse beam (with 
a period of t—152 as). In the middle row, the superposition 
of spherical harmonics with widths given by the numerical 
simulation is repre sent ed; they show perfect agreement with 
the selection rules (4j5l. (b) Same as upper row in (a) but for 
a left-circularly polarized pulse beam. 



that 52% of the quantum electron state is ionized when 
the beam is linearly polarized, and 31% when the beam 
is circularly polarized. Note that Eq. ^ yields a to- 
tal photoionization probability that depends non-lincarly 
on the light polarization. We can express the electron 
quantum state at each time as \ip(t)) = a\ipi) + \5ip(t)), 
where Sip is the excited state part. The excited state 
function can be decomposed in an unbound spherical ba- 
sis dip(r) = Y^l M u L,M( r )Yif (fit f)- O ur initial elec- 
tron state has full spherical symmetry, belonging to the 
spherical harmonic Y®. However, after the interaction 
with the pulse, the electron state excites different spher- 
ical harmonics. In Fig. [3] the projection of the excited 
state onto the xy-plane is depicted during the first cycle 
of the pulse beam, being a superposition of spherical har- 
monics obeying the selection rules Q. Also, depending 
on the input polarization, the electron evolution varies 
noticeably. 

By resorting to numerical approaches, we could accu- 
rately evaluate the projections of the excited states onto 
the spherical harmonics Yff and extract the correspond- 
ing probabilities Pl,m = / dr r 2 \uL t M{r)\ 2 ■ Using this 
numerical method, the widths of the spherical harmonic 
superpositions at different times have been derived, show- 
ing excellent agreement with the electron state evolution, 
as represented in Fig. [3]ja). Moreover, we have ana- 
lyzed the lowest spherical harmonic content of the final 
excited electron state (see Fig. [1]) in three scenarios: (i) 
with a Gaussian pulse beam (in the transverse spatial 
approximation, it could be considered as a plane wave) 
linearly polarized in the x-direction (the electron state 
is ionized about 30%); (ii) the case for a pulse spatially 
modulated by a Laguerre-Gaussian mode, linearly polar- 
ized in the cc-direction; (iii) when the Laguerre-Gaussian 




FIG. 4: Final Quantum Electron State and their 
Spherical Harmonic Spectrum. Projection of the ex- 
cited state (J dz\5ip(x,y, z)\ 2 ) onto the plane xy after the 
interaction with the pulse beam for: (a) a Gaussian mode, 
(b) a Laguerre-Gaussian mode linearly polarized and (c) a 
Laguerre-Gaussian mode circularly polarized, (d) Spectrum 
of spherical harmonics for the three different cases (a), (b) 
and (c). 



pulse is left-circularly polarized. The main remark is that 
the spherical harmonics Yf 1 and Y{~ are most efficiently 
excited by the plane-wave-like pulse, in striking contrast 
with the Laguerre-Gaussian scenario, where no such ex- 
citation exists. If the Laguerre-Gaussian pulse is linearly 
polarized, then, Y® , Y% and F 4 4 are the most occupied 
states, whereas if it is circularly polarized, F ° and Y§ are 
the most relevant. There is a small contribution from Y®, 
as the electron is less ionized. We emphasize that these 
results are in accordance with the derived selection rules 
@ and (5). 

It is interesting to examine the total angular momen- 
tum transferred to the quantum electron state. First, we 
calculate the mean value of the electron OAM, (L) = 
(r x (p — q A)), during its evolution. Figure [5] shows the 
time evolution of the orbital angular momentum along 
the ^-direction in the same case of Fig. [3ja), while a 
depopulation of the fundamental state occurs. The elec- 
tron starts in the ground state, with zero OAM. As the 
pulse begins to interact with the electron, the OAM of 
the latter in the z-direction oscillates, but notice that at 
the end of the pulse, the electron quantum state gains 
a finite amount of OAM: 1.53 au (1.53 ft). There is no 
OAM contribution in other directions. We expect, as 
the ground state has no OAM in the absence of a field, 
that the electron excited states belong to unbound states 
bearing OAM. On the other hand, when the pulse is left- 
circularly polarized, as the case of Fig. ^b), the OAM in 
the z-direction is negligible, as it is expected. Regarding 
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FIG. 5: Temporal Evolution of the Electron Angular 
Momentum. Solid red curve represents the OAM of the 
electron in the z-direction during interaction, parameters as 
in Fig. |3|a). Notice that the electron begins with null OAM 
and after the pulse beam it gains up to 1.53 au (1.53 ft). Dash 
gray curve represents the population of the ground state. As 
time increases, it is ionized up to 0.52 of the initial population. 



the excited state position components, mean values are 
zero except for the z-component, where a small shift of 
10~ 2 au is present, caused by a non-vanishing magnetic 
field at the origin. 



IV. DISCUSSION 



to achieve high-harmonic generation [2H US] with OAM, 
the experimental challenges to extend this phenomena to 
more complex systems (such as Rydberg atoms) and the 
possibility to use OAM for nuclear quantum optics ap- 
plications [26, 27 . For instance, by exploiting the trans- 
verse profile, one could address the Ml transition at 3.5 
eV of Th-229 [2S]. 



V. APPENDIX 

We derive in this section the selection rules for the 
interaction of a light beam possessing orbital angular 
momentum with matter. From now on, it is conve- 
nient to take the quantization axis in the beam prop- 
agation direction, the position vector can thus be writ- 
ten as r = (x,y,z) — r (sin 9 cos (j>, sin 9 sine/), cos 9) and 
p = r sin 9. Hence, the interaction element given by equa- 
tion (13]) is proportional to 
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A J dr r 3 f^y^l u L f ,M f {r)u Lz , Mi (r) x 
[a (sin0) l<l+1 cos0e- i(ut -^- x) +c. c]y%' + 



We present the first work addressing the photoion- 
ization process induced by beams carrying OAM. We 
have found novel selection rules Q and ^ for a pulse 
beam characterized by a topological charge I. In addi- 
tion, other interesting effects have been revealed. If the 
electron excited state after interaction is let to evolve, 
it goes away from the origin. In contrast, by intro- 
ducing a pulse beam with N cyc = 14, we notice that 
the excited electron state remains confined, within a ra- 
dius of 10 au, due to the ponderomotive force induced 
by the Laguerre-Gaussian profile. We also simulate the 
case where the atom is displaced 2 au from the origin in 
the x and y-directions, with N cyc — 14. In both cases, 
the atom is ionized much faster since the electric field is 
more intense now, and the excited electron state remains 
trapped. From evaluations of the mean values of the 
position component, we have observed an electron mo- 
tion around the vortex. Varying the polarization of the 
beam, in particular for left-circular polarization, an ion- 
ized state with a ring structure (AM — 0) is predicted, 
see Figs. [3^b) anQc). Furthermore, modifying the ini- 
tial phase, different ionized state structures can be gen- 
erated. Thus, by tunning the phase and the polarization, 
one expects a manipulation of the ionized state. There 
still remain many open questions, such as the feasibility 



(6) 

Equation ^ can be divided into two parts; one depends 
exclusively on the radial contribution and the other one 
on the angular contribution. If we expand the angular 
contribution of Eq. ^ , 

(sin^l+icos^e-^ 4 -^*) = e ~' ( ^~ x) x 
{(sm0)W +1 e^ +1 ^ + (sin^l+V^ 1 ^} , 

(sin^l+^in^e-^*-^-*) = * 
{(sin#)l £ l +1 e^ +1 ^ - (sin^l+V^ 1 )*} , (7) 

we can notice basically two terms that are repeated 
( sm 0)KI+V(^+i)</> an d (sin 0)1*1+ V^" 1 )*. These two 
terms can be decomposed in spherical harmonics func- 
tions, getting then the selection rules. In order to obtain 
the spherical harmonic decomposition, we must consider 
two cases, when I > and I < 0. Let us begin with case 
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I > 0, where the angular terms read as 
(sin^l+V( |+1 ^ = (_i)^ 2 ( £+1 )(£ + 1)! x 



4tt 



Y t 



y(2£ + 3)! 1+1 
(sin^l+V^" 1 ^ = {-l)hV- l \t - 1)! 



(8) 
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(2£-l)! 



where we have used sin 2 9 = A\fh{Y§ - Y§ /v / 5)/3 and 
# = (-1) V(2^ + 1)! sin* e i ^/2^!V47T. The first term 
Q, in integral pi), gives rise to the straightforward se- 
lection rule: |A£| < * + 1, AL + t + 1 is even and 
AM = £+1. In equation (|9| there are two contributions, 
the first one yields also straightforwardly the selection 
rule: | AL\ < £ - 1, AL + £ - 1 is even and AM =1-1, 
but the second one is a bit subtle. First of all, we need to 
decompose the product Y^ZiY^ m t° spherical harmon- 
ics. Using the following formula; 



Y m i y™ 12 

h h 



^ 2 g J (2l 1 + l)(2l 2 + l) x 

i=|/l— (2! m= — l 



4tt(2Z + 1) 

(li,/a;0,0|I,0)<li,Za;m 1 ,m 2 |i ) m> If , (10) 

where (Zi, Z2; rni, "12^3, fia) are the corresponding 
Clebsch-Gordan coefficients, we know that it is possible 



to write Yf~}Y? = aY 



bY 



cY e _ 1 + dY l _ 2 



-1 x 2 — Lt J e+i ~i~ u x t 
e YfZ§- We can take cZ = e = 0, as M < L must be satis- 
fied in a spherical harmonic function. And 6 = always, 
due to (£ - 1, 2; 0, 0\£, 0) = 0. Concretely, 



' {21- 


1)5 


Att{2£ 


+ 3) 


1 {21- 


D5 



{£- 1,2; 0,0|*+1,0) x 
(£-l,2;i-l,0|* + !,*-!) , 



(*-l,2;0,l 



1,0) x 



y 4tt(2£- 1) 

(£— 1,2;€ — 1,0|£ — — 1> . 

Hence, the second contribution from Eq. ([9| to the selec- 
tion rules can be summarized as: |AL| < £+1, AL + ^+1 
is even and AM = I — 1. 

Now, let us consider the case when I < 0. Proceed- 
ing in an analogous way than before, we can now write 
expressions (|8| and ^ as: 



(sin^l+V^ 1 ^^^- 1 ) 



47T 



(2|*|-1)! 



(sine)W+ 1 e^- 1 ^ = 2^l+ 1 )(|f| + l)!x 

/ 47T 



(2|*| + 3) 



T Kl+i 
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where we have used sin 2 9 
Y[ = J{2\£\ + l)\sm w 6 1 
term (]12[). in integral (6j), yields: |AL| < 




4 v / 5(r o ° - F 2 °/V5)/3 and 
/2l^|^|!\/47r. The second 
+ 1, 

1 is even and AM = 1—1. On the other 
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hand, in Eq. (11) there are two contributions, the first 
one giving rise to: |AL| < \£\ - 1, AL + \£\ - 1 is 
even and AM = £ + 1. In the second one, the product 
Y^P^Y^ must be decomposed into spherical harmonics 
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using formula (10 1. As before, we can write Y^^ 1 ^ 



■c'Y, g+1 
L I \t\-x 



d'Y t+1 



e' Y^} 3 , where 



dl = e' = owing to |M| < L, which must be satisfied 
in any spherical harmonic function. And b' = always, 
due to (|*| - 1, 2; 0, 0| |£| , 0) = 0. Therefore, 
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-1,2;£ + 1,0||£|-1,£ + 1) . 



implying that: |AL| < \£\ + 1, AL + \£\ + 1 is even and 
AM = * + 1. 

Generalizing all the last calculations, the selection 
rules derived from the matrix element ^ for beams car- 
rying any unit of orbital angular momentum are summa- 
rized in Eq. Q. 

At variance with the plane waves £ = {\AL\ = 1), we 
can expect larger exchange of angular momentum. Of 
course, playing with the beam polarization (a and /?), 
as we can note in equation Q, we can modify the selec- 
tion rules. For example, for a right-circular polarization 
(a = 1 and /3 = i), the only surviving terms are given 
by expressions (|8|) and (11), restricting AM = £ + 1. 



Analogously for a left-circular polarization {a = 1 and 
P = — i) the only surviving terms are given by expres- 
sions j9| and (12), restricting AM = £ — 1. Selection 
rules pT, in photon terms, can be thought as the absorp- 
tion of a photon carrying a total angular momentum in 
the propagation direction m = £+s, where s indicates the 
polarization part (spin momentum, for right-circular po- 
larization s = 1 and for left-circular polarization s = —1). 
We would like to remark that these selection rules are 
exclusive to the transverse profile. Moreover, the second 
interaction hamiltonian TH.ii, in the case of plane waves, 
is just a constant term, yielding a ponderomotive force. 
In our case, it is quite different. Analogously to equation 
([3]), we can write 



(*f\Hnm = ± n (* f \A 2 {r,t)\* i ), (13) 

in order to extract the selection rules. Again, considering 
the vector potential (1), in which the dipolar and trans- 
verse spatial approximation has been taken into account, 
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we can write 



yielding the selection rules of Eq. |5]) . 



(#/|«zz|*i) oc 



(tf/lfAopl'l e-^-W-x) + c . c.l |#.). (14) 



Eq. (14 1, if we expand the quadratic potential, two dis- 



tinguishable terms appear: 
A o/9 l z l e -K^-^-x) + c c l 2 = 

A 2 r 2 l'l [(a 2 + /? 2 )(sin0) 2 l £ l e -«(^-<*-x) + c . c . 

+(A P W ) 2 , 



where the second term does not depend on time, it acts 
as a well-potential. On the other hand, the first term 
gives rise to new selection rules as 



(sin0) 2 IV 2 ^ = 2 2e (2£)\. 



4tt 



{4£ + l)\ 



Y 



21 



VI. ACKNOWLEDGMENTS 



We acknowledge support by the Spanish Ministry of 
Education and Science under contracts FIS2005-01369, 
FIS2008-02425, FIS2006-04151, FIS2007-29091-E, and 
Consolider projects SAUUL and QOIT, CSD2007-00013, 
CSD2006-00019, and the Catalan, Junta de Castilla y 
Leon, and Junta de Castilla-La Mancha Governments 
under contracts SGR2005-00358, SA146A08, and PCI08- 
0093. We also acknowledge R. Corbalan for fruitful dis- 
cussions. 



L. Allen, M. V. Beijersbergen, R. J. C. Spreeuw, and J. 
P. Woerdman, "Orbital angular momentum of light and [12 
the transformation of Laguerre-Gaussian laser modes," 
Phys. Rev. A 45, 8185 (1992). 

A. T. O'Neil, I. MacVicar, L. Allen, and M. J. Pad- [13 
gett, "Intrinsic and Extrinsic Nature of the Orbital An- 
gular Momentum of a Light Beam," Phys. Rev. Lett. 88, [14 
053601 (2002). 

S. Franke-Arnold, L. Allen, and M. Padgett, "Advances [15 
in optical angular momentum," Laser & Photonics Re- 
view 2, 299 (2008). [16 
M. Bhattacharya and P. Meystre, "Using a Laguerre- 
Gaussian Beam to Trap and Cool the Rotational Motion 
of a Mirror," Phys. Rev. Lett. 99, 153603 (2007). 
M. F. Andersen, C. Ryu, P. Clade, V. Natarajan, A. [17 
Vaziri, K. Helmerson, and W. D. Phillips, "Quantized 
Rotation of Atoms from Photons with Orbital Angular [18 
Momentum," Phys. Rev. Lett. 97, 170406 (2006). 
A. Mair, A. Vaziri, G. Weihs, and A. Zeilinger, "En- [19 
tanglement of the orbital angular momentum states of 
photons," Nature 412, 313 (2001). 

R. Inoue, N. Kanai, T. Yonehara, Y. Miyamoto, M. [20 
Koashi, and M. Kozuma, "Entanglement of orbital angu- 
lar momentum states between an ensemble of cold atoms 
and a photon," Phys. Rev. A 74, 053809 (2006). 
G. F. Calvo, A. Picon, and A. Bramon, "Measuring two- [21 
photon orbital angular momentum entanglement," Phys. 
Rev. A 75, 012319 (2007). 

M. van Veenendaal and I. McNulty, "Prediction of Strong [22 
Dichroism Induced by X Rays Carrying Orbital Momen- 
tum," Phys. Rev. Lett. 98, 157401 (2007). 

I. G. Mariyenko, J. Strohaber, and C.J. G.J. Uiterwaal, [23 
"Creation of optical vortices in femtosecond pulses," Opt. 
Express 13, 7599 (2005). 

I. Sola et al., "High power vortex generation with volume [24 
phase holograms and non-linear experiments in gases," 



Appl. Phys. B 91, 115-118 (2008). 

A. Einstein, "Concerning an Heuristic Point of View To- 
ward the Emission and Transformation of Light," Ann. 
Phys. 17, 132 (1905). 

P. B. Corkum and F. Krausz, "Attosecond science," Na- 
ture Phys. 3, 381 (2007). 

C. I. Blaga et al., "Strong-field photoionization revis- 
ited," Nature Phys. 5, 335 (2009). 

G. A. Mourou, "Optics in the relativistic regime," Rev. 
Mod. Phys. 78, 309 (2006). 

J. R. Vazquez de Aldana and L. Roso, "Magnetic field 
effects in strong field ionization of single-electron atoms: 
Three-dimensional numerical simulations," Laser and 
Part. Beams 20 185-193 (2002). 

S. M. Barnett, "Optical Angular-Momentum Flux," J. 

Opt. B: Quantum Semiclass. Opt. 4 S7 (2002). 

S. J. Van Enk, "Selection rules and centre-of-mass motion 

of ultracold atoms," Quantum Opt. 6, 445 (1994). 

R. Jauregui, "Rotational effects of twisted light on atoms 

beyond the paraxial approximation," Phys. Rev. A 70, 

033415 (2004). 

A. Alexandrescu, D. Cojoc, and E. Di Fabrizio, 
"Mechanism of Angular Momentum Exchange between 
Molecules and Laguerre-Gaussian Beams," Phys. Rev. 
Lett. 96, 243001 (2006). 

G. F. Calvo, A. Picon, and R. Zambrini, "Measuring the 
Complete Transverse Spatial Mode Spectrum of a Wave 
Field," Phys. Rev. Lett. 100, 173902 (2008). 
Shawn. A. Hilbert et al., "Temporal lenses for attosecond 
and femtosecond electron pulses," Proc. Natl. Acad. Sci. 
106, 10558 (2009). 

T. Brabec and F. Krausz, "Intense few-cycle laser fields: 
Frontiers of nonlinear optics," Rev. Mod. Phys. 72, 545 
(2000). 

S. Patchkovskii, Z. Zhao, T. Brabec, and D. M. Vil- 
leneuve, "High Harmonic Generation and Molecular Or- 



8 



bital Tomography in Multielectron Systems: Beyond the 
Single Active Electron Approximation," Phys. Rev. Lett. 
97, 123003 (2006). 

[25] J. Courtial, K. Dholakia, L. Allen, and M. J. Pad- 
gett, "Second-harmonic generation and the conservation 
of orbital angular momentum with high-order Laguerre- 
Gaussian modes," Phys. Rev. A 56, 4193 (1997). 

[26] K. W. D. Ledingham, P. McKenna, and R. P. Sing- 
hal, "Applications for Nuclear Phenomena Generated by 



Ultra-Intense Lasers," Science 300, 1107 (2003). 
[27] T. J. Biirvenich, J. Evers, and C. H. Keitel, "Nuclear 

Quantum Optics with X-Ray Laser Pulses," Phys. Rev. 

Lett. 96, 142501 (2006). 
[28] E. Peik and Chr. Tamm, "Nuclear laser spectroscopy of 

the 3.5 eV transition in Th-229," Europhys. Lett. 61, 181 

(2003). 



